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Abstract 

This paper shows that there is a correspondence between quasi-exactly solv- 
able models in quantum mechanics and sets of orthogonal polynomials {P n }- 
The quantum-mechanical wave function is the generating function for the 
P n (E), which are polynomials in the energy E. The condition of quasi-exact 
solvability is reflected in the vanishing of the norm of all polynomials whose 
index n exceeds a critical value J. The zeros of the critical polynomial Pj(E) 
are the quasi-exact energy eigenvalues of the system. 
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In quantum mechanics there exist potentials for which it is possible to find a finite 
portion of the energy spectrum and associated eigenf unctions exactly and in closed form. 
These systems are said to be quasi-exactly solvable. [Ij-Q In such systems the potential 
depends on a parameter J; for positive integer values of J one can find J eigenvalues and 
eigenfunctions exactly. The usual approach to the analysis of quasi-exactly solvable systems 
is an algebraic one in which the Hamiltonian is expressed as a nonlinear combination of 
generators of a Lie algebra, not belonging to the center of the corresponding enveloping 
algebra. [§J This technique is a modification of the dynamical symmetry approach to exactly 
solvable quantum-mechanical systems, in which one can find by algebraic means the entire 
spectrum in closed form. 

In this paper we propose an alternative approach to quasi-exact solvability. We show 
that the solution ip to the Schrodinger equation for a quasi-exactly solvable model, 

Hip = Eif), (1) 

is the generating function for a set of polynomials {P n (E)} in the energy variable E. These 
polynomials satisfy a three-term recursion relation and therefore form an orthogonal set 
with respect to some weight function w(E). For positive integer values of the parameter J, 
corresponding to quasi-exact solvability, we find that the norm of P n {E) vanishes for n > J. 
Moreover, all polynomials P n (E) beyond a critical polynomial Pj(E), factor into a product 
of two polynomials, one of which is Pj(E): 

P n+J (E) = Pj(E)Q n (E) (n>0). (2) 

The zeros of the critical polynomial Pj{E) are precisely the quasi-exact energy eigenvalues 
of the quantum-mechanical model. 

We illustrate these features of quasi-exactly solvable models with the following infinite 
class of Hamiltonians first discussed by A. Turbiner |IJ 
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Here, s is an arbitrary parameter. When s lies between ~ and |, there is an attractive 
centrifugal term; for s outside this range the centrifugal term is repulsive. When s = j or s = 
|, the centrifugal core term disappears leaving a nonsingular sextic oscillator Hamiltonian 

H = — - (4s + 47 - 2)x 2 + x 6 . (4) 
ax 2 

When the parameter J in Eq. (H) is a nonnegative integer, the corresponding Schrodinger 
equation has J exact, closed-form solutions for any value of s. 

We seek a solution ip(x) to the Schrodinger equation for H in Eq. of the form 

W .«p(-^)^e(-I)"^Lj»-. (5) 

Observe that when s = ~ this solution becomes an even-parity wave function of the oscillator 
Hamiltonian (|J); when s = |, ^(x) becomes an odd-parity wave function of in (f|). 

Demanding that i/j(x) in Eq. @ obey the Schrodinger equation (0) leads to the following 
recursion relation for the expansion coefficients P n (E): 

P n (E) = EP n _i(E) + 16(n-l)(n-J-l)(n + 2s-2)P n _ 2 (E) (n > 2), (6) 

subject to the initial conditions 

P Q (E) = 1 and P X {E) = E. (7) 

From these initial conditions the recursion relation (|6]) generates a set of monic [|7| polyno- 
mials, the next four of which are 

P 2 {E) = E 2 + (32 - 32J>, 

P 3 (E) =E 3 + [(160 - 96 J)s - 32 J + 64}E, 

P 4 (E) =E 4 + [(448 - 192J)s - 128 J + 352]£ 2 + (3072 J 2 - 12288 J + 9216)s(s + 1), 
P 5 (E) =E 5 + [(960 - 320J> - 320 J + 1120]£ 3 + [(15360J 2 - 81920J + 91136)s 2 

+ (25600 J 2 - 141312J + 164864)s + 6144J 2 - 36864J + 49152]E. (8) 

These polynomials have a number of noteworthy properties. First, for all values of the 
parameters s and J they form an orthogonal set. This follows from the fact that they 
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are generated by a second-order (three-term) recursion relation. || The appearance of a 
three-term recursion relation is a consequence of the form of the potential in Eq. (|3|). For 
example, the corresponding recursion relation for an x A anharmonic oscillator potential, 
whose Hamiltonian is not quasi-exactly solvable, is a higher-order recurrence relation. The 
harmonic oscillator system leads to a two-term recursion relation; this system is exactly 
solvable rather than quasi-exactly solvable. 

Second, from the expansion we can see that the wave function ip(x, E) is the gener- 
ating function for the polynomials P n (E). 

The third and most significant property of the polynomials P n {E) is that, when the 
parameter J takes positive integer values, the polynomials exhibit the factorization property 
in Eq. (|2|). This factorization occurs because the third term in the recursion relation (||) 
vanishes when n = J + 1, so that all subsequent polynomials have the common factor 
Pj(E). This factorization property holds for all values of the parameter s. Furthermore, 
this factorization leads to the result that the zeros of the critical polynomial Pj{E) are just 
the quasi-exact energy eigenvalues. This is true because the expansion in (|JJ) truncates when 
E is a zero of Pj(E); when this series truncates the wave function ip(x) is automatically 
normalizable. 

To illustrate this factorization we list in factored form the first six polynomials P n {E) 
for the case J = 3: 

Po(E) = 1, 

Pi(E) = E, 

P 2 (E) = E 2 — 64s, 

P 3 (E) = E 3 - (128s + 32)£, 

P 4 (E) = [E 3 - (128s + 32)E]E, 

P 5 (E) = [E 3 - (128s + 32)E}(E 2 + 128s + 192). (9) 
Observe that Ps{E) is a common factor of P n (E) for n > 3. The zeros of P^{E) are 

E = 0, E = ± v / 128s + 32, (10) 
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which are the three exact energy eigenvalues for the quasi-exactly solvable Hamiltonian ([3D 
when J = 3. The corresponding exact eigenf unctions are obtained by evaluating ip(x) in 
Eq. (pi) at these values of E: 



i>o(x) = ex P ( — x 4 } 1 - 

V K ' F \ 4 / r(2s) V 2s + 1 

= exp f-^ 4 ) 



1 4 \ x 25 - 1 / 2 / V128S + 32 2 x 4N 

f(2sfV 8i x + 2~ S/ 

1 A x 25 - 1 / 2 / Vl28s + 32 2 , x v 
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Note that the energy levels may be ordered by the number of nodes of the corresponding 
wave function. 

A fourth property of the polynomials P n (E) concerns their norms. The norm (squared) 
7„ of P n (E) is defined as an integral: 



J dEw(E)[P n (E)] 2 . (12) 



It is possible to determine the norms of an orthogonal set of polynomials directly from the 
recursion relation; it is not necessary to know explicitly the weight function w(E) with 
respect to which the polynomials are orthogonal. || The procedure is simply to multiply 
the recursion relation @ by w(E)E n ~ 2 and to integrate with respect to E. Using the fact 
that P n (E) is orthogonal to E k , k < n, we obtain a simple, two-term recursion relation for 

In- 

7„ = 16n(J — n)(2s + n — 1)7„_ 1 . (13) 

The solution to this equation with 7 = 1 is 

16 n n\T(J)T(2s + n) 
ln ~ T(J~n)T(2s) • ( } 

This equation reveals that the space of orthogonal polynomials arising from a quasi- 
exactly solvable model is associated with a nonpositive definite norm. In particular, we can 



see from Eq. ([TJ|) that 7„ vanishes for n > J if J is a positive integer. The appearance 
of a vanishing norm coincides with the factorization mentioned above and is an alternative 
characterization of quasi-exact solvability. 

It is interesting that while the polynomials P n+ j(E) for n > have vanishing norm when 
J is a positive integer, the quotient polynomials Q n (E) in Eq. (g) form a new orthogonal 
set of polynomials for each value of J. 

Having determined the norms 7„ of the polynomials P n {E) it is natural to evaluate the 
integral of the square of the generating function (wave function) with respect to the weight 
function w(E): 

G(x) = J dEw(E)[^(x } E)] 2 , (15) 

where ip(x,E) is given in Eq. fl5|). Using the orthogonality of the polynomials P n (E), we 
can express G(x) as a confluent hypergeometric function: 

r(j) / 1 ^ » x 4 ^ 4 *- 1 . 

G{X) = TWf XP \T2 X ) S n!T(n + 2s)T(J-ny (16) 

When J is a positive integer, this sum truncates and we find that G(x) can be expressed as 
a linear combination of the squares of the J quasi-exact eigenfunctions of the Hamiltonian 
H in Eq. (§). For example, when J = 3, we have 

/ i \ i i _|_ 



s 2s(2s + l) 

= r(2s) {Wi [Mx)]2 + 4TTi [Mx)]2 + 4lTT [ ^ (x)]2 ) ' (17) 

where ipo(x) and ip±(x) are taken from Eq. (|ll]). We emphasize that this result is highly 
nontrivial. Expressing G(x) as a linear combination of the squares of the eigenfunctions 
requires that one solve an overdetermined system of 2J — 1 equations for J expansion 
coefficients. 

Let us now investigate the properties of the weight function w(E). From the polynomials 
P n (E) we can calculate the moments of w(E). Let a n represent the 2nth moment of w(E): 




dEw(E)E 



(18) 



(Because the polynomials have parity symmetry we know that the odd moments vanish.) 
We are free to normalize w(E) so that its zeroth moment is unity: 



The remaining moments can then be determined algebraically: 
ai = 32(J- l)s, 

a 2 = 32 2 (J - l)s(3Js - 5s + J - 2), 

a 3 = 32 3 ( J - l)s(15 J 2 s 2 - 60Js 2 + 61s 2 + 15 J 2 s - 67 Js + 74s + 4J 2 - 19 J + 22), 
a 4 = 32 4 (J - l)s(105J 3 s 3 - 735J 2 s 3 + 1743Js 3 - 1385s 3 + 210J 3 s 2 - 1596J 2 s 2 

+4038Js 2 - 3372s 2 + 147 J 3 s - 1179J 2 s + 3114Js - 2688s + 34J 3 - 282J 2 

+765J- 674), 

a 5 = 32 5 (J - l)s(945J 4 s 4 - 10080J 3 s 4 + 40950J 2 s 4 - 74400Js 4 + 50521s 4 + 3150J 4 s 3 
-35910J 3 s 3 + 153990J 2 s 3 - 292154Js 3 + 205228s 3 + 4095 J 4 s 2 - 48960 J 3 s 2 
+218337JV - 427524Js 2 + 307860s 2 + 2370J 4 s - 29306J 3 s + 134373J 2 s 
-269085 Js + 197206s + 496 J 4 - 6272 J 3 + 29292J 2 - 59531J + 44134). (20) 

These moments have some interesting mathematical properties. For example, all the mo- 
ments a n , n > 1, have a factor of (J — l)s. Furthermore, in the residual factor the coefficient 
of (Js) n_1 is (2n — 1)!! and the coefficient of s n_1 is the nth Euler number E n . [|l(j] 

The outstanding property of the moments a n concerns their rapid rate of growth. This 
rate of growth can be determined using the fact that there is a simple relationship between 
the moments a n and the coefficients of P n _2{E) in the recursion relation (|6|). Specifically, 
the Taylor series 



a = 1. 



(19) 



CO 



/(*) = E 



n 



(21) 



n=0 



whose coefficients are the moments in Eq. (0), is equivalent to a continued fraction 
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f(z) = 1/(1 - hz/(l - b 2 z/(l - b 3 z/(l - . . .)))), (22) 



whose coefficients are b n . 11 Since b n is a cubic polynomial in n we deduce that the moments 



a n grow like (3n)\. p~T 



It is unusual to find orthogonal polynomials whose weight functions have moments that 
grow so rapidly. The classical orthogonal polynomials, such as the Hermite polynomials, 
typically have moments that grow like n\. This is also true of discrete versions of the clas- 
sical orthogonal polynomials, such as the Hahn polynomials. [|12| The Euler and Bernoulli 
polynomials are distinctive |13] in that their moments grow like (2n)!. However, the poly- 



nomials P n {E) associated with quasi-exact solvability are of an entirely new type due to the 
rapid rate of growth of their moments. Carleman's condition states that when the moments 
grow faster than (2n)!, the moment problem is not guaranteed to have a unique solution. [15| 



Almost certainly, the weight function w(E) is not unique! This nonuniqueness corresponds 
to a kind of gauge invariance that underlies these quasi-exactly solvable systems. Indeed one 
may conjecture that the nonuniqueness of the weight function is related to the Lie algebraic 
symmetry of quasi-exact solvability. 
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